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Abstract — We face the following dilemma for designing low- 
density parity-checls codes (LDPC) for quantum error cor- 
rection. 1) The row weights of parity-checli should be large: 
The minimum distances are bounded above by the minimum 
row weights of parity-check matrices of constituent classical 
codes. Small minimum distance tends to result in poor decoding 
performance at the error-floor region. 2) The row weights of 
parity-check matrices should not be large: The sum-product 
decoding performance at the water-fall region is degraded as 
the row weight increases. 

Recently, Kudekar et al. showed spatially-coupled (SC) LDPC 
codes exhibit capacity-achieving performance for classical chan- 
nels. SC LDPC codes have both large row weight and capacity- 
achieving error-floor and water-fall performance. In this paper, 
we design SC LDPC-CSS (Calderbank, Shor and Steane) codes 
for quantum error correction over the depolarizing channels. 

Index Terms — spatial coupling, LDPC code, iterative decoding, 
CSS codes, quantum error-correcting codes 

I. Introduction 

In 1963, Gallager invented low-density parity-check 
(LDPC) codes III, which is defined as a kernel of a sparse 
parity-check matrix. Due to the sparseness of the parity- 
check matrix, LDPC codes are efficiently decoded by the 
sum-product (SP) algorithm, messages of SP decoding can be 
statistically evaluated. Since the optimized LDPC codes can 
approach very close to the Shannon limit |2|, error-correcting 
code theorists are attracted by LDPC codes. 

By the discovery of CSS (Calderbank, Shor and Steane) 
codes 1 3 1, |4| and stabilizer codes [51, the notion of parity- 
check measurement is introduced as a generalized notion of 
parity-check matrix. From a point of view of this general- 
ization, quantum LDPC codes are naturally defined via low- 
density parity-check measurements |j6] by Postol. In particular, 
a parity-check measurement for a CSS code is characterized 
by a pair of parity-check matrices. If both of these parity- 
check matrices of the pair are sparse parity-check matrix, the 
related CSS code is called a quantum LDPC (CSS) code. 

For classical case, a random construction method for con- 
structing a sparse parity-check matrix generate a high error- 
correcting performance LDPC code. On the other hand, it is 



almost impossible to apply the same method for quantum case, 
since the pair of parity-check matrices have to satisfy the 
following constraint: the product of one of the pair and the 
transposed other is subjected to be a zero-matrix. Therefore 
one of the research interests for quantum LDPC codes is stated 
as: "Achieve the constraint and sparseness simultaneously. 
Additionally, valuable structure for classical LDPC codes is 
also achieved." 

MacKay et al. proposed the bicycle codes [7| and Cayley 
graph based CSS codes ITOl . These codes are known as 
self-dual containing LDPC codes. In fSl, two of the authors 
proposed a construction method of CSS code pair that has 
quasi-cyclic ( QC) parity-check matrices with arbitrary regular 
even row weight dr > 4 and column weight d\ such that 
rfi /2 > di > 2. In [91, the authors generalized this construction 
method to codes over non-prime field of characteristic number 
2. To the best of the authors' knowledge, these codes (T\, 
ifTOl . ifTTj hold the highest error-correcting performance among 
efficiently decodable quantum LDPC codes so far. 

Spatially-coupled (SC) LDPC codes are classical capacity- 
achieving codes based on the construction of convolutional 
LDPC codes. Felstrom and Zigangirov llT2l introduced a con- 
struction method of {d\, (ir)-regular convolutional LDPC codes 
from (di, (ir)-regular block LDPC codes. Surprisingly, the 
LDPC convolutional exhibited better decoding performance 
than the underlying block LDPC codes under a fair comparison 
with respect to the code length. Note that in this paper, 
convolutional LDPC codes are the LDPC codes defined by 
band sparse parity-check matrices. 

Kudekar et al. named this phenomenon as "threshold satu- 
ration" and proved rigorously for the binary erasure channels 
fT3\. Further, Kudekar et al. 1 14 1 observed empirical evidence 
which support the threshold saturation occurs also for the 
binary-input memoryless symmetric-output (BMS) channels. 
Another remarkable advantage of SC LDPC codes is encoder 
universality. In other words, keeping di/dr, in the limit of 
large dj, L and w, the coupled ensemble {d\,di-, L,w) ifTSl 
achieves universally the capacity of the BMS channels under 



SP decoding. Conventional capacity-achieving codes such as 
polar codes (TS] and irregular LDPC codes fZ] do not support 
such universality. According to the channel, polar codes and 
need selection of frozen bits ifTTl and irregular LDPC codes 
need optimization of degree distributions. 

In this paper, we propose a construction method for spatially 
coupled quantum LDPC codes for quantum error correction. 
By the proposed method, we obtain a pair of parity-check 
matrices such that these matrices are each orthogonal and each 
matrix has a diagonal band structure. 

II. Preliminaries 

Let us recall definitions of related fundamental notion of 
quantum LDPC codes and classical SC LDPC codes. 

A. CSS codes 

A CSS code Q, which is the main interest of this paper, 
is a class of quantum error codes. The code is a complex 
vector space. The vector space Q is characterized by a pair 
of classical binary linear codes C and D whose parity-check 
matrices He and Ho satisfy HcHj, — 0. 

A CSS code Q associated with {Hc,H]j) is defined as a 
complex linear combination of the following vectors: 



E 



d!) for c e C, 



where is the dual code of as a classical code and the ba- 
sis for quantum states is assumed as a standard computational 
basis. 

B. LDPC-CSS codes 

When the parity-check matrices He and Hd are sparse, the 
CSS code associated with {He, Ho) are called an LDPC- 
CSS code. As it is written above, the structure of the code 
space Q is completely characterized by two linear codes C 
and D associated with {Hq^Ho). The aim of this paper is 
to construct these codes as LDPC-CSS codes, particularly SC 
LDPC codes. 

One of issues for designing [He, He,) is the row weights 
of these matrices. LDPC-CSS codes are efficiently decoded 
by using SP syndrome decoding [7J. Its success probability 
for decoding is highly affected by the row weight of these 
parity-check matrices {Hc,Hd)- 

1) The row weights of each of He and Hd should be 
large: The minimum distances of C and D are less than 
or equal to the minimum row weights of Hjj and He, 
respectively. Small minimum distance tends to result in 
poor decoding performance at the error-floor region. 

2) The row weights of each of He and Ho should not 
be large: The SP decoding performance at the water- 
fall region are degraded if row weight of Hd and He 
increase. 

In summery, we are facing the dilemma above for designing 
{He, He,)- The proposed idea to overcome the dilemma is to 
introduce a class of classical "spatially coupled LDPC codes" 
to quantum LDPC codes. 



C. Sparse Band Model Matrices and Spatially Coupled LDPC 
Codes 

Here, we construct classical SC LDPC codes from sparse 
band model matrices T-L. Let d\,dt and ric be positive 
integers. The integer d\ shall be the column weight of Ti.. 
Let Tig be a positive integer such that rig divides d\. Briefly 
speaking on the structure of H, submatrices of size d\ x dt 
appear ric times on the diagonal line. Thus the column weight 
of H is di. Define d^ := dtd\/ns- The row weight of H 
is dr in the middle and dt in the top and bottom. Define 
M := di + {ric — l)ns and N := ricdt- The size of H is 
M X N. 

We give a formal construction for the sparse band model 
matrix. Define a matrix H = (hjj) over {0, *} of size {d\ + 

(ric — l)"-s) X ncdt by putting hj i := * if 6ns + 1 < i < 
{b + l)ns or hj i = 0, where * is a formal symbol and b := 

La-i)Mj. 

The following is an example of a sparse band model matrix 
H for parameters d\ — A,dt — 2, ric = 12 and Ug = 1. 

**00 

******** 
******** 
******** 
wwwwyy******** 

11= Q Q Q Q Q Q Q Q s i * * * * * * 

""******** 
******** 
******** 
'******** 
****** 
* * * * 
00**, 

The sparse band model matrix H contains nc(= 12)-subblocks 
of size di{= 4) x dt{= 2). 

Next, we define a {di, , ric) SC LDPC code via the model 
matrix Ti. Let P be a positive integer. Replace each "*" of "H 
with a binary permutation matrix of size P x P and each 
of H with a zero matrix of size P x P. The replaced matrix 
is called spatially couple LDPC matrix and its kernel space is 
called a classical spatially coupled LDPC code. For fixed dr 
and d\, the size PM x PN of such a sparse band matrix scales 
as M = Ndi/d,. + di - 1 increases. The value R:^ l-M/N 
is called a design rate. The design rate converges to 1 — d\/dr 
as ric tends to infinity. The SC LDPC code ensemble lfT3l 
was shown to exhibit SP decoding performance which is very 
close to the MAP decoding performance of (di, dt)-regular 
LDPC code ensemble ifTSll in the limit of large code length 
and coupling number ric- For a sparse band matrix H, define 
a matrix H' = {hj i) over {0,*} by putting h'^ i 




hr- 



where r is the number of rows of H. We also call H' a sparse 
band model matrix. 

SC LDPC codes enable to have large minimum row weight 
of parity-check matrix and simultaneously exhibit excellent 
SP decoding performance. Hence, in this paper, we study the 
design of SC LDPC-CSS codes. It is expected that such codes 
overcome the dilemma on the row weight. 

D. Quasi-Cyclic LDPC Codes 

In the previous works, the decoding performance of a SC 
LDPC code is analyzed by the technique based on ensemble 



analysis. In other words, the instance of SC LDPC code shall 
be constructed from randomly chosen permutation matrices 
for a sparse band model matrix. On the other hand, it is not 
expected that the requirement for {Hc,Hd), i.e. HcHj) ~ 0, 
holds on parity-check matrices of SC LDPC codes by using 
random choice. Therefore, we slightly relax the constraint of 
spatially-coupled LDPC codes. We employ circulant permuta- 
tion matrices instead of permutation matrices. 

Let us define a matrix /(I) over binary field of size P x P 
by putting: 



10 

10 

■■. 

1 

1 



G{0,1} 



PxP 



In other words, /(I) is a circulant permutation matrix which 
is a single right shifted matrix for the identity matrix. For an 
integer x, define a P x P matrix I{x) :— /(l)^. I{0) is the 
identity matrix. For a symbol oo, define a P x P matrix /(oo) 
as a zero matrix. 

For a matrix (cj,;)o<j<di,o<;«i, over {0, 1, . . . , P- 1, oo}, 
let H he a Pd\ x Pdi- binary parity-check matrix defined as 
follows: 

H := {I{cjj))o<-i<di,o<e<dr- 



We refer to such a matrix as (d\ , , P) quasi-cyclic (QC) 
permutation matrix, or simply a QC permutation matrix. 

The associated linear code is called a QC LDPC code. 

III. Spatially Coupled LDPC-CSS Codes based on 
Quasi-Cyclic Permutation Matrices 

Let {Hc,Hd) be an LDPC-CSS code. If He and Hd are 
SC LDPC matrices, we call the LDPC-CSS code a spatially 
coupled LDPC-CSS code. It is known that small cycles in 
the Tanner graph degrade the SP decoding performance. In 
this section, we construct SC LDPC matrix pair [He, Ho) 
which satisfies the following three conditions. 

1) HcHl = 

2) Both of the Tanner graphs of He and Hd are free of 
cycles of length 4, 

3) He and Hd are QC-permutation matrices, 

A. Conventional QC-LDPC CSS codes 

We first review the construction developed in fSl. The 
authors proposed [8J the following method for constructing a 
QC parity-check matrix pair {Hc,Hd)- In the original paper 
JS), the construction method is flexible about the row size of 
the matrices, i.e.. He and Hd can have different row sizes. 
For simplicity, in this paper, we focus our attention on the 
case that He and Hd have the same row weight d^ and have 
the same column weight d\. 

Let [cj^i) and (dfc.z) be matrices over {0, 1, . . . , P — 1} U 
{oo} of size d\ x d^. Let He (resp. Hd) be a QC-permutation 



matrix associated with (cj,;) (resp. {dk,i))- From ||8l, it is 
proved that He Hi = 0, if #{£ G [0, P - 1] | c^- ^ - dui = p 
mod P} is even for all j E [Q,d\ — 1] and for all p e [0, P — 1] . 
The Tanner graph of He is free of cycles of length 4, if 



#{c,,£ - ck.i eZ*p\£e[0, d, - 1]} = d. 



(1) 



for aU p e [0,P - 1], and for afl j,k e [0,di - 1] such 
that j ^ k. We give the following theorem which is slightly 
generalized version of Theorem 6.1 in |8l. 

Theorem III.l. Let P be a positive integer with P > 2. 
Define Zp := {z £ Zp | 3a G ZpjZa = 1}. For a,T E Zp, 
define ord(cr) :— min{m > | ct™ = 1}, and define (r)^ = 
{r, Tcr, ...,Tcr°'-^('^)-l}. 
Let di, be integers and ti, T2 be in Zp such that 

d\ > 2,d, > 4, 

dr/2 = ord(f7), (2) 
1 < di < ord((T), 
ord(a) ^ #Z*p, 

1 - e Zp /or all 1 < j < ord(cr), 

T2^(<7),,. (3) 

Let He and Hd be two {di,d,-, P)-QC binary matrices such 
that 

He = {I{cjj))o<j<duQ<e<d^, 
Hd = {I{dj,e))o<j<di,o<£<dr, 



where 



Tia 

T2 



-j+e 



o<e< dj2 
0--^'+^ d,/2<e < d„ 



-T2(7- 



Q<1< d,/2 
d,/2 <£< dy. 



Then, it holds that 1) HeHj) = 0, 2) both of the Tanner 
graphs of He and Hd are free of cycles of length 4. 

See for the construction proof. 

If Ti — 1, the statement is the same as the one in 
From Theorem lIII.il we obtain two diPx d^P binary matrices 
He and Hd such that HeHj, — and the Tanner graphs of 
He and Hd are free of cycles of size 4. In this paper, we 
refer these codes as conventional QC-LDPC CSS codes. The 
following is an example: 



Example III.l. With parameters d\ ^ 3, d^- = 6, P = 7, a = 2 

and Ti = 1,T2 =3, from Theorem \IIL1\ we obtain a JP x 
di-P binary matrix pair {He^Ho) such that HeHj) — as 
follows. 
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Fig. 1. Illustration of parity-check matrices Hq and //jj of a spatially-coupled LDPC-CSS code (C,D). 



B. Spatially coupled CSS codes 

Let Hc{ti,T2) and Hd{ti,T2) denote a pair of parity- 
check matrices associated with ti and T2 in Theorem IIII. ll 
Let Uc pairs of matrices Hc{t\^\t2^) and iJ^) (r}*-' , t^*') 
for i — 0,...,nc — 1. Let rig be a positive number such 
that ns\d\. We construct sparse band QC matrix pair He and 
Hd as depicted in Fig. [T] It holds that HcHJ) = since 

HciT['\4'^)HD{T^'\4'^V = for z = 0,...,nc - 1. 
Therefore the orthogonaHty HcHjj = does not depend 
on the choice of t^^^ , T2*^ . He and ffc have regular-column 
weight d\ and slightly irregular row weight. The row weight is 
typically dtd\/ns around the center of He and Ho and takes 
the minimum weight dt at the boundaries. 

On the other hand, the condition "there are no cycles of 
size 4 in the Tanner graph" depends on the choice of them, 
since the Hc{ti\t2^) and Hc{t[^ \'''2 "*) share a common 
row if \i — i'\ < di/ug. By a similar argument in the proof of 
Theorem lIII.il we have the following theorem: 

Theorem m.2. // (r^)^ n {t^'')^ = % for any h,h' e {1,2} 
and i, i' G {0, 1, . . . , ric — 1} such that \i — i'\ < di/ug, then 
HqHJ) — and there are no cycles of size 4 in the the Tanner 
graphs of sparse band QC matrices He and Hd. 

Example III.2. In Fig. |2] we give an example of He and Hu 
of SC LDPC-CSS codes with d\ = 3, dt = 6, P = 31, = 

6,^3 = l,(rr,rf)) - (16,4),(r(^\Tf') = (8,12), 
(rf ,rf)) = (6,l),(rf ,rf)) = (3, 11), (r^^), )) = 
(17,2),(rf ,rf)) = (6,4). 

Remark IILl. The Tanner graphs for He and Hjj are iso- 
morphic. This statement is proved by combinatorial argument 
although we omit the proof in this paper 

IV. Numerical Results 

We assume the transmission takes place over depolarizing 
channels Q Section V] with depolarizing probability 2/ni/3, 
where /m can be viewed as the marginal probability for X 
and Z errors. Note that the channel is the normal depolarizing 
channel. We assume the decoder knows the depolarizing 
probability 3/ni/2. 



Figure |3] compares the proposed SC LDPC-CSS codes and 
conventional QC-LDPC CSS codes. The quantum coding rate 
of these codes is 0.50 and 0.49, respectively. Such rate-loss is 
due to the coupled construction. The rate-loss can be reduced 
by increasing ric. The parity-check matrices of these codes 
have the same column weight 10. The minimum row weight 
40 and 20, respectively. For both proposed and conventional 
CSS codes, two component classical codes C and D are 
isomorphic. This isomorphism implies that C and D have the 
same average decoding performance. Hence, we restrict our 
attention only to the performance of C. 

Due to the large minimum row weight of parity-check 
matrix, all the codes exhibit no error-floors at the bit error 
rate down to 10~^. No undetected errors were observed for 
both codes. Hence, it is expected the minimum distances of 
both codes are sufficiently large and tight up to the upper 
bounds 40 and 20, respectively. As posed in the dilemma, it is 
observed that conventional QC-LDPC CSS codes have only a 
small coding gain by increasing code length from n = 100840 
to n = 412840. The proposed code outperforms the conven- 
tional code with much shorter code length n — 101000 and 
n = 181000. This excellent performance outweighs the rate- 
loss. 

V. Conclusion 

In this paper, we proposed spatially-coupled LDPC-CSS 
codes. Numerical experiments show that the proposed codes 
have both deep error-floors and excellent water-fall perfor- 
mance over the depolarizing channels. In other words, thresh- 
old saturation is also observed for the LDPC-CSS codes. 
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Fig. 2. Parity-check matrices Hq, Ho of the proposed SC LDPC-CSS codes defined by dj = 3, dt = 6, P = 31, ric = 6, ns = 1. /(■) is abbreviated for 
each entry. 
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Fig. 3. Comparison of decoding en'or rate of the proposed SC LDPC-CSS 
code (red) with {d\ = 10, dt = 20, rtc = 50, ris = 5) and conventional QC- 
LDPC CSS codes (blue) with (d\ = 10, dr = 40). The quantum coding rate 
is Rq = 0.49 and Rq = 0.50, respectively. The code length is n qubits. 
The block (dark) and bit (light) en'or rate of the constituent codes C and D 
of the proposed CSS code pair (C, D) over the depolarizing channel with 
marginal flip probability /m of X and Z eiTors. Due to the isomorphism, the 
error rates of C and D are the same. 
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VI. Appendix 
Proof for SC Construction (Theorem UII.li : 1) Note that 

{c,,, - c,,,, e Zp I < £ < |} = (ri(a-^- - a-^'))^, 

{c,,, - c,v,, eZp\^<e<d,} = {t2{<j-^ - <y-'')),, 

where e [0, c?i — 1] such that j j' . Therefore 

#{Cj- f - CfJ G Zp I < ^ < y } - d,/2 

and 

#{cj,e - c/.f eZp\^<e<d,.}= d,/2. 
In general, 

{x)„ n {x')a = ^ {xy)^ n {x'y)„ ^ 
for any x,x' ,y £ Zp. By (O, 

(Ti)^n(r2)^ = 0. 

It implies that 

{n{a-^ - a-^'))^ n (T2(a-J- - a"-'"')), = 
^ (ri), n {t2), ^ 0. 



Hence 



#{c,,^-c,',f eZp \ o<e<d,} 

-C,',, eZp |0<^< y} 

+ #{Cj,£ - Cy,f e Zp I y < £ < 

dj2 + dj2 = d,.. 



By ([T]i, there are no cycles of size 4 in the the Tanner graph 
of He- We can show for Hjj in a similar way. 
2) By direct calculation, 

{cjj - dk.i e Zp I < ^ < y } 

= {cj/ - dk,i e Zp I y < ^ < d,.}. 

Therefore HcHj, = 0, if #{£ G [0, F - 1] | c^- £ - dfe,^ = p 
mod P} is even for all j G [0, d\ — 1] and for all p G [0, P— 1]. 
It means HcHj, =0. ■ 



